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, ■ Abstract 



A number of internal inconsistencies, misleading conclusions and lack of 
completeness in the analyses of the small x deep inelastic data in recent papers 
of HERA physicists, particularly in a very recent one by the HI collaboration 
are pointed out. It is also shown that, when an analysis without prejudice is 
^ ' carried over, the situation is very different from the one claimed there in what 

•/^ . respects checks of theoretical expectations. 

^ . 

■ The so-called HI Collaboration has recently produced a preprint concerning deep inelastic data at 

I small X, that we will henceforth call HERA2. The work constitutes an improvement of the former HERA 

analysisl^J of the structure function F2, to be denoted by HERAl, extending it to wider x, ranges and 
greatly diminishing the errors. 

There is not much to be argued about the experimental obtaining of the data. The machine and the 
pi il detectors worked wonderfully. The quality of the data is excellent. It is therefore more of a pity that the 

^ ■ physicists in the collaboration have not been more careful with the subsequent theoretical analysis which 

shows features that, to put it as mildly as possible, are difficult to justify. Because the second paper HERA2 
(ref.2) is so much superior in the quality of the data, I will mostly discuss it here. Regarding it, I will address 
k>( \ myself to a number of separate (but related) questions which I will now ennumerate: 1) Parametrizations 

. of the data. 2) Behaviour of structure functions as a; — > 0. 3) Theoretical analysis with so-called "Double 

d ' scaling limit". 4) Connection with various theoretical models. This is the central part of the present note, 

which will then be finished with some comments. 

l.-The theory of strong interactions is QCD. On this practically everybody is agreed, so, when presenting 
parametrizations of data one should try to, at least, not violate grossly standard QCD features. In this, I 
think, it is no excuse that the parametrization is merely "a phenomenological ansatz" , as is described the 
parametrization of p. 20 of HERA2, Eq. (8). As should be by now widely known, the structure function 
F2 may be split in a singlet and a nonsinglet piece, F2 — Fs + F^s- For small x, Fj^ig ~ B^s^'^'^- 
matter what, the fact that both Fg and Ff^s ^^rs proportional to cross sections implies that Bjqg > 0: 
"phenomenological ansatz" or not. A measure of the quality of the fit is obtained by noting that the 
corresponding coefficient in (8), denoted by e in the table underneath this Eq. (8) in HERA2 is negative. 
Really. 

But there is more. It may perhaps be argued that Eq.(8) was not meant to have anything to do with a 
QCD description: at least it should be compatible with the rest of the paper HERA2. In the same page 20, 
the authors state that the wide range of covered by their experiment allows them to study the behaviour 
at small x, and its variation with Q^. So they write 

F2(x,Q2) ^ ^-^ (1) 



> 
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and adscribc this prediction to Do Riijula ct alJ'"*! Wc will come to this last point later; for now what interests 
us is that the HERA2 collaboration find, and report in Table 4, values of A increasing as* 

A = 0.24 = 12 GeV^) to A = 0.50 (Q^ = 800 GeV^). (2) 

Actually, a careful look at Table 4 (p. 23) shows that, at least up to = 350 GeV^, A may be considered 
constant modulo what may easily be interpreted as statistical fluctuations, with a value of A around A ^ 0.32. 
Again, we will return to this later; for now what interests me is the blatant contradiction between Eq.(2) 
and the parametrization of Eq.(8) where, for x — > 0, one has F2 — cx'^ with d = —0.188, widely off the 
bounds given in Eq. (2) . How the authors of HERA2 may claim that their data are fitted, at the same time 
by a;^^ '^^ and x^''^'^'^ baffles me. But there is more. 

2. -It is extremely important, for judging the quality of the theoretical analysis of HERA2, to return to 
the matter of behaviour at small x. The authors of HERA2 fit F2 at small x; and they are so certain of their 
results that they display them not only in the text, but in Table 4 (to which I have already referred). The 
smallest value of A they find is for = 2.5 GeV'^, where they get A = 0.19; for = 1.5 GeV^ the value is 
0.21. In fact, an analysis at = (Compton scattering) gives W a A around 0.26. 

Then the authors of IIERA2 state that "the rate of growth of F2 is expected to increase" , and promptly 
signal to ref.3 as the place where such growth was predicted. So one opens ref. 3, and there one reads 
the prediction (p. 1651, line 9ff): "The rate of growth increases with increasing Q^. It is always weaker 
than a power'^ of x, although stronger than a log. (Italics mine). So, the HERA people who both in their 
phenomenological parametrization [Eq. (8) of HERA2] and in their Table 4 find a power, conclude that they 
check the prediction of ref. 3 -a surprising conclusion. 

3. -To clarify matters a bit more, I would like to bring attention to the derivation of the results of ref.3. 

Specifically, in p. 1650 these authors consider the moments of the structure fimctions. They realize [their Eq. 
(2)] that the moments are given in terms of the Wilson expansion which, after a slight change of their rather 
old fashioned notation to the one prevalent nowadays, implies the relation 

) ^ — 77^ an ^ - — an. (3) 
The a„ may easily be identified with expectation values of operators. 



an ~ {p\On\p) 

that need not be specified here. The D{n) are related to the anomalous dimensions. De Riijula et al. 

recognize that the behaviour of F2 for small x is linked to the singularities (in n, considered as a continuous 
variable) of the moments; the leading behaviour being given by the rightmost singularity. This is clear since 
one has ^ 

Mn((3')= / dxF2{x,Q^)x^-\ 
Now, if one assumes that at a certain Qq one has a Pomeron-type behaviour, 

F2{x,Ql) ^ Co, Co = constant, (4) 

cc— >0 

then one has, at that Qq, a singularity for n = 1 that one can identify with that of the anomalous dimension, 
D(n), which indeed is singular at n = 1. 

It should be emphasized that here one has a dichotomy. If at a given Qq one has a behaviour like 
a constant, then the structure function, by virtue of the analysis of De Riijula et al.I^l will grow "always 
weaker than a power" of x, in their own words. Contrarywise, if at any Q"^, F2{x, Q^) grows like a power, 

* We only report the values from = 12 GeV^; below this value the contamination of nonsinglet is very 
strong, and even rough analyses shoud incorporate this fact to get any credibility. 
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the conditions of the theorem cannot apply. Actnahy, the situation is really very simple. As discovered long 
agoI^J, one has two possibilities. Either the singularities of the matrix elements a„ lie to the left of those of 
D{n), i.e., to the left of n = 1, or to the right of it. If the first, one has the situation envisaged in ref.3. 
If, however, the dominating singularity of a„ lies to the right of n = 1, at a certain n = 1 + A, A > 0, then 
one has a behaviour like a power, for all Q^, F2{x,Q'^) ~ The coefficient of proportionality, a function 
of Q^, may in fact be calculated. One may describe this result by stating that in the low x regime, the 
dependence on x and of structure functions "factorizes" . 

The paper of De Riijula et al. is not very explicit. To find a more detailed calculation, we will refer to 
the review of Ball and Forte'^l, apparently known to the HERA2 collaboration since they quote it. There, in 
p. 9, the authors explicitely state that, if F2 ^ x~^, the chain of arguments which lead to a "double scaling 
form", namely their Eq.(2.31), breaks down, and a power-like behaviour becomes dominant [their Eq.(2.34), 
when one undoes their changes of variable, reads exactly F2{x,Q'^) ~ a;~^]. Curiously enough, this power 
behaviour is exactly what is found by both HERA collaborations. 

We can be more specific. Assume, with De Rujula et al, and Ball and Forte, that, contrary to the 
findings of the HERA2 group, one had, at a given Qq, a Pomeron-like behaviour as that given by Eq.(4) in 
the present paper. Then, as correctly stated in refs.3, 5, one obtains a behaviour dominated by the (known) 
singularities of D{n). A simple calculation then gives, for all ^ Qq, x — > 0, the very explicit formula 



F2{x,Q^) 



Co 



33 -'2nj- 



5767r2|loga;|log[a,(g2)/a,(Q2 



exp - 



' 144| logx| 
(33 -2n/) 



log 



Q.(Oo) 



asm 



(5) 



This equation is essentially identical to the corresponding one in ref.5, p. 9; the behaviour it implies has been 
described as "double asymptotic scaling" . Roughly (but we will come to more precise statements below) Eq. 
(5) states that 



^2(2;, Q^) ~ exp v^l logx|/i. 



(6) 



On the other hand, and in the same notation, the findings of the HERA2 fits, in exponential form, imply 

F2(a;,Q2)~exp{|loga;|/2}, (7) 

where the fi may depend weakly on Q"^. Undaunted by the evident incompatibility of (6) and (7), the 
authors in HERA2 happily assert at the same time the validity of (7), in their Table 4, and Eq. (8), and the 
validity of (6) in p. 26: "Thus double asymptotic scaling is a dominant feature in this region". 

4. -How can these contradictory statements be reconciled? They cannot. The physicists in HERA2 have 
failed, in their theoretical analysis of the data, to do what an honest experimentalist should do, viz. to test 
theoretical assumptions without bias. 

In fact there existed at the time of the first and (of course) the second HERA analyses three theoretical 
predictions (that I know of). Prom a standard moments analysis one can either have the "double asymptotic 
scaling", refs.3, 5: 



F2{x,Q^) Co 



33 -2n/ 



5767r2|loga;|logK(Q2)/a,(Q2)] 



exp ^ 



' 144|loga;| 
(33 - 2n/) 



log 



2m-Dii^0.5 



(8) 



This depends on the singularity of D{n) being the rightmost one. On the other hand, if it is the singularity 
of a„ that lies to the right, one gets the factorization behaviour given in ref.6: 



F2ix,Q^) Bs[as{Q^)]-<'+x-^ + BNs[as{Q^)]-''''x°-'; 



(9) 



I have added the nonsinglet contribution -Bjvs'[as(Q^)] Dn ~ 0.512, to both expressions for future 

ease of reference. In (9), d+ is the largest eigenvalue of D, easily calculated; an explicit formula for it may 
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be found in rcfs.4 or 6. It is to be stressed that, in the present state of the art, one cannot decide on the 
basis of perturbative QCD alone* on which is vaUd, (8) or (9): it should be for experimentalists to tell us. 

There is yet a third theoretical estimate [''1. It is very difficult to compare with the other two, as it is 
obtained using totally different methods to get the prediction 

F2{x,Q') ivo = Ca = 3. (10) 

TT 

Let me finish this point by telling what one gets if trying to fit (8), (9) or (10) to the HERA data. This 
test will be done in two steps. I will choose only data with x < 10^^, to be sure that one is in the "small 
a;" region and, in the first fits, I will also take 100 GeV^ > Qq > 10 GeV^ to avoid problems with varying 
number of excited flavours, large NLO corrections, and to have, for the test of (8), ^ Qo, whatever the 
last may be. This region contains a total of 48 high-quality points. In the second step, the set of points will 
be enlarged to up to 350 GeV^. We will then have 63 experimental points, distributed over a wide range 
in Q^. Thus, the analysis will necessitate inclussion of NLO corrections. To the approximation needed, these 
may be implemented by simply taking as{Q'^) to two loops, and allowing A to vary freely, which is what I 
will do here. 

The results are now summarized, starting with the first situation. (10) runs contrary to the trend of 
the data, and produces a very large chi-squared, x^/d.o.f. = 926/(48 — 2). What is more, the fit rejects 
the introduction of a NS component, a clear indication of the fact that the behaviour (10) is certainly not 
attained at the HERA energies (the fit to tlw^ higher data does not fare much better, either). Thus we will 
say no more about (10): the interested reader may find further discussion, and a possible re-interpretation 
of (10) in ref.4. 

For the fits with Eqs.(8), (9), we fix n/ = 4, take to one loop, 

, 2^ 127r 

as[Q )- (33_2„/)iogQ2/^' 

and choose A = 0.2 GeV so that as{m'^) = 0.32. The results would not change much qualitatively if allowing, 
e.g., yl to be a free parameter (some such results may be found in ref.4, and in the present note below). 

The fit with the formula of De Riijula et al., not taking into account the NS piece is quite bad, but the 
results improve when including a NS contribution, as is done in Eq. (8). They are summarized in Table Al. 

Qi B^s Qf> xVd-o.f. 
7.1 X 10-3 2.0 0.10 GeV^ 51.3/(48 - 3) 

Table Al.- Fit with "double scaling", Eq.(8). A = 0.2 (fixed) 

This should, of course, be compared with the fit one gets (Table Bl) if using the equation (9), since (9) 
and (8) stem from mutually incompatible assumptions. If we neglected the NS contribution in (9), we would 
also get a poor chi-squared; when including the NS piece this improves to about one unit by d.o.f., see Table 
Bl. 

Bs Bns a xVd.o.f. 
1.04x10-2 1.21 0.40 49.8/(48-3) 

Table Bl.- Fit with "factorization", Eq.(9). A = 0.2 (fixed) 

The results reported in Table Bl are more satisfactory than those given in Table Al: not so much the 
chi-squared by d.o.f., only marginally better, but because of the following other features. Firstly, the value 
of Bns in Table Bl is closer to what one obtains in the analysis of structure functions like uW^ in neutrino 
scattering, which are pure nonsinglet and which suggest B^s ^ 0.6. Another problem with the results 
reported in Table Al is that they imply that the structure function at Qq, essentially proportional to the 
Compton scattering cross section, should behave as a constant with the energy (proportional to 1/x), which 
it does not: as shown in refs.4. 8, 9, it still grows like a power of 1/x. 

* For deep inelastic -yp scattering. For 7*7 the situation is different, and will be discussed below. 
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The situation is even more clear if we use the second set of HERA2 data with .t > 10 ^, and all 
between 12 and 350 GeV^. Here one has to take into account NLO corrections, which, as explained, may be 
well approximated by just using the two- loop formula for as{Q'^), leaving yl as a free parameter. 

For the "double scaling" hypothesis, Eq.(8), the resulting chi-squared/d.o.f. deteriorates. Numerically, 
I have conducted a fit with (8) and to two loops and find the results of Table A2. 

A Co B^s Ql xVd-o.f. 

0.40 GeV 4.9x10-3 2.0 0.38 GeV^ 69.9/(63-4) 

Table A2.- Fit with "double scaling", Eq.(8); A taken as a free parameter, n/ = 4. 

That (8) does not give a very brilliant fit may also be seen in Fig. 10 of HERA2 (p. 25), particularly in the 
lower part. In spite of the use of a log log plot, on which anything looks like a straight line, eyeball inspection 
of the scattering of the experimental points around the theoretical curve indicates a chi-squared/d.o.f. clearly 
larger than unity (the authors of HERA2 seem to be allergic to numbers, and do not give any figure for this 
or several other chi-squarcd of their fits). On top of it, the two problems mentioned in connection with the 
former fit, non-Pomeron behaviour at Qq = 0.38 GeV^ and deviation of Bats from its expected value persist. 

For Eq.(9), the details of the analysis may be found in ref.4; the chi-squared per d.o.f. actually improves 
to less than one unit, and B^s decreases to the value of 0.8, quite close to the neutrino result of 0.6. This 
is summarized in Table B2. 





A Bo 


Bns 


A xVd.o.f. 




0.10 GeV 1.0x10- 


-3 0.77 


0.38 50.4/(63-4) 


Table B2.- 1 


Fit with "factorization'' 


. Eq.(9): 


.1 taken as a fvvc parameter. /) /■ = 4. 



The x^/d.o.f. is now substantially better in the fit with (9) then in the fit with (8), 50.4 vs. 69.9; but 

other features also improve. Thus the value of B^s in Table B2 is now comfortably close to that expected 
from neutrino deep inelastic scattering, 0.77 compared to 0.6. 



5.- Comments. One can of course still not not draw definite conclusions, but there is a clear trend: 
the data of HERA, both the old HERAl and the new HERA2, support clearly a "factorization" behaviour 
at small x, as predicted on the basis of the analysis of lower energy data long ago by Lopez and the present 
authorial. As shown in the present note the "double scaling" description, which the paper HERA2 claims to 
provide evidence for* ("double asymptotic scaling is a dominant feature") actually fails to give a good fit to 
the data, besides presenting extra problems. It is really sad that the HERAl and HERA2 collaborations have 
not taken the trouble to present a fair analysis of experiment, not such a difficult task as I have been shown 
here. Why they ignore the possibility of a power- like behaviour as that given in Eq.(9) is really difficult 
to understand; apart from the original papers, it is given in the textbook of ref.lO considered (by e.g. the 
Particle Data Group) as one of the standard references in QCD. The situation is even more outrageous since 
many of the physicists involved in HERAl and HERA2 have also been involved in the PETRA analyses of 
deep inelastic j*j scattering for which Wittenl^^l has derived (and experiment has checked) a behaviour, 
at small x, exactly like that in (9), with a value of A, A ~ 0.5 constant up to a slight dependence in the 
number of fiavours excited. This was indeed one of the theoretical reasons given in the original papers^^l for 
preferring (9) to the "double scaling" behaviour. 



* In all fairness it should be noted that the authors of HERA2 admit the existence of some deviations of 
the data from the predictions of double asymptotic scaling, e.g., in the last line of p. 24. 
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